Dielectric and thermal relaxation in the energy landscape 
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We derive an energy landscape interpretation of dielectric relaxation times in undercooled liquids, 
comparing it to the traditional Debye and Gemant-DiMarzio-Bishop pictures. The interaction be- 
tween different local structural rearrangements in the energy landscape explains qualitatively the 
recently observed splitting of the flow process into an initial and a final stage. The initial mechan- 
ical relaxation stage is attributed to hopping processes, the final thermal or structural relaxation 
stage to the decay of the local double-well potentials. The energy landscape concept provides an 
explanation for the equality of thermal and dielectric relaxation times. The equality itself is once 
more demonstrated on the basis of literature data for salol. 

PACS numbers: 64.70.Pf, 77.22.Gm 



I. INTRODUCTION 

Broadband dielectric spectroscopy [l| is the most ver- 
satile method to study the flow process in molecular liq- 
uids. However, the quantitative nature of the relation 
between dielectric signal (i.e. molecular reorientation) 
and shear flow is not yet clear. The classical Debye pic- 
ture and its extension to viscoelasticity [2, LJ considers 
the molecule as a small sphere with a hydrodynamic ra- 
dius th immersed in the viscoelastic liquid. It predicts a 
slow dielectric decay, about a factor of fifty slower than 
the mechanical shear stress decay. The measured dielec- 
tric decay is on the average a factor of ten faster than 
this Debye prediction A thorough quantitative anal- 
ysis of dielectric and shear data in seven glass formers 
IE showed a general qualitative agreement with the 
Gemant-DiMarzio-Bishop extension 0, 0] of the Debye 
scheme, but a rather poor quantitative fit. Very recently, 
an alternative to the Debye model and its extensions was 
proposed Q, which gave a much better fit for glycerol 
y], |g and propylene carbonate Q data. 

The proposal was based on the growing evidence for 
the identity of thermal and dielectric relaxation functions 
[^ IhIIT^.Ii3| (dielectric hole burning experiments 0,0] 
even suggest that dynamically distinct domains in the 
liquid are associated with a time constant characterizing 
both the dielectric and the thermal behavior). 

The present paper begins with a reminder of the text- 
book introduction to dielectrics |l(J in terms of the de- 
cay of the dielectric polarization after a switch-off of the 
electric field. In the Debye scheme, the decay of the 
dielectric polarization is due to the rotational diffusion 
of the molecules with a Debye-Stokes-Einstein diffusion 
constant. But one can easily generalize the formalism to 
other decay mechanisms. In particular, we consider the 
recent proposal pj of an initial and a retarded part of 
the flow process. We show that one must expect such a 



division in the energy landscape picture. The thermal re- 
laxation sees both the initial and the retarded final part 
of the process. The equality of thermal and dielectric 
relaxation times is understandable in terms of an energy 
landscape argument. The equality itself is checked for 
salol, using literature data. 



II. THEORETICAL BASIS 



In SI units, the electrostatic equations read 



D Q = e o e(0)E Q = e E Q + P a 



(1) 



where eo is the vacuum permittivity, Eq is the electric 
field, Do is the displacement field and Po is the polar- 
ization per unit volume. Eq, Dq and Po are vectors, so 
the static dielectric constant e(0) is a second rank tensor. 
Here, however, we limit ourselves to isotropic liquids or 
glasses, so e(0) reduces to a scalar. Similarly, we neglect 
the conductivity contribution, assuming that it can be 
simply subtracted. 

The electrostatic equations are generalized to 
frequency-dependent equations p~(j 



D(lo) = e e(w)P(w) = e E(uj) + P{uj) 



(2) 



with frequency-dependent fields and dielectric constant, 
respectively. Thus one has 
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The next step introduces a useful simplification for , 
the dielectric constant at infinite frequency. Since one is 
mainly interested in the relaxation in the frequency range 
below 100 GHz, one considers all processes above this fre- 
quency as immediate processes. This includes the elec- 
tronic polarizability (time scale 10 -15 seconds) as well as 
the vibrational contributions (time scale 10 -12 seconds), 
so 
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where n is the refractive index and Ae V ib is the contri- 
bution of the vibrations (molecular librations) to the di- 
electric constant. Since this contribution is nonzero, one 
expects > n 2 in a reasonable fit (this is in fact a main 
problem of the Gemant-DiMarzio-Bishop extension 0, 
of the Debye scheme 0). 

Consider an electric field Eq staying constant from t = 
—00 until t = 0. At time zero, the field is switched off 
(Fig. 1). The Fourier transform of this field reads 



-Eq. 



(5) 



For such a field, the polarization has a constant value 



P Q = (e(0) - l)e E 



(6) 



for times t < 0. At t = 0, the polarization drops instan- 
taneously to Pq — Poo , where 



Poo = (foe - l)eo^o 



At t > 0, P(t) decays with 



P(i) = (P - Px=)$(t), 



(7) 



(8) 



where $(i) is a function which begins with 1 at time zero 
and drops to zero at infinite time. 

If one compares this dielectric polarization decay with 
the mechanical shear stress decay, one finds a fundamen- 
tal difference. After switching off the electric field, the 
polarization decays without any externally applied force. 
In the mechanical case, one applies a small external shear 
strain at time zero, keeps the strain constant and one 
observes the decay of the shear stress. The mechanical 
shear stress decay is described in linear response by a 
time-dependent shear modulus G{t). The dielectric case 
is different. In order to describe it in terms of a time- 
dependent dielectric constant, one needs to go to the 



time 



opposite case of Fig. 1, switching the electric field on 
at time zero and watching the rise of e(t) from to e 
with the function 1 — $(t). As a consequence, if one wants 
to compare decay with decay in order to see which quan- 
tity decays faster, one needs to compare the mechanical 
modulus with the dielectric susceptibility and not with 
any dielectric modulus, at variance with a recent proposal 
|5J. We come back to this point after the treatment of 
the Gemant-DiMarzio-Bishop relation 0, [j| . 
The Fourier transform P(ui) of P(t) is given by 



P(w) = Po + (Po - Poo)P[$(t > 0)], 



(9) 



where F[$(t > 0)] denotes the Fourier transform of the 
decay function Q(t). 

Inserting eqs. J5j and (JSJ into eq. |JSJ, one gets 

e^) = l + -(^ + ^-^ l u J FMt >())}). (10) 

The relation between $(cj) and <&(i) is the same as the 
one between G(u>) and G(t) 

P 00 

$(w) ee -iwF[$(t > 0)] = u) I $(t)(smut+icosujt)dt. 

Jo 

$(w) is a complex function which is 1 for infinite fre- 
quency and decreases to zero as the frequency goes to 
zero. 

With this definition, and inserting eqs. JfJJl and JJJ into 
eq. Q1U[). we arrive at the final result 
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(12) 



The Debye decay mechanism is the rotational diffusion 
of the molecules, with the diffusion constant given by the 
Dcbye-Stokes-Einstein equation 



D t 



k B T 4 2 
7, = ^r H D rot , (13) 



where D trans is the translational diffusion constant of the 
molecule and D ro t is its rotational diffusion constant, th 
is the hydrodynamic radius of the molecule (note that the 
hydrodynamic Debye-Stokes-Einstein equation is in prin- 
ciple derived for larger objects than a single molecule). 
For continuous rotational diffusion, the relaxation time 
for the Legendre polynomials is 



1 



T~L,rot 



L{L + l)D ro t 



(14) 



where L is the order of the Legendre polynomial. For the 
dielectric signal, L = 1, one obtains the Debye relaxation 
time 



FIG. 1: Polarization decay after switch-off of the electric field 
at time zero. 



(15) 
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The decay function for the dielectric polarization is a 
single exponential (Debye process) 



which leads to the well-known Debye equation 



6(0) 



1 + iUT D 



(16) 



(17) 



To extend the Debye equation to viscoelasticity 0, |j| , 
one replaces the static viscosity r\ by a frequency- 
dependent function r\{io\ which is in turn related to the 
frequency-dependent shear modulus G(iS) 



iu>rj(uj) = G(uj) = Goog{u>). 



(18) 



Here Goo is the infinite frequency shear modulus and g(u>) 
is a normalized function like $>(io), going from zero to 1 
with increasing frequency. 

The ratio between viscosity and infinite frequency 
shear modulus defines the Maxwell time 



tm 



Gnc 



(19) 



the characteristic shear relaxation time. 

If we replace the viscosity rj in the td of the right hand 
side of the Debye equation by the expression of eq. (|18fl , 
we get the Gemant-DiMarzio-Bishop expression 



1 



1 



1 + (4 7 rG 00 r| f /fcflT) ff (a;) 1 + c r g{oj) ' 



(20) 



which contains the dimensionless ratio cy between the 
Debye relaxation time tjj and the Maxwell time tm- 

The Gemant-DiMarzio-Bishop expression, eq. (|20|l . 
does not approach zero for infinite frequency as the De- 
bye expression in eq. I|17|) . because the molecule is still 
able to turn in a balance between the torques exerted 
by the electric field and the finite restoring force of the 
sheared elastic medium, respectively. In fact, one has in 
this case 



Ae,„;h — 



e(0) 



1 + cy 



(21) 



Thus Ae v ib is already taken into account, and the 
Gemant-DiMarzio-Bishop equation |2j, |3j takes the form 



1 



:(0) - n 2 1 + c r g(uj) ' 



(22) 



With a hydrodynamic radius of 0.2 nm, a Goo of 2 GPa 
and a temperature of 200 K one calculates c r = 50, 
so one expects a small vibrational component and nearly 
two decades difference between td and the Maxwell time. 
The first of these expectations is in fact found, but the 
second fails by an average factor of ten 0, 0, S US • 

Note that the failure is only in the factor and not in 
the temperature dependence. Unlike the translational 
diffusion, which decouples from the viscosity below some 



critical temperature [l6l Il7l Hsf . the rotational relaxation 
time follows essentially the temperature dependence of 
the viscosity over the whole temperature range. 

On the basis of the Gemant-DiMarzio-Bishop equation 
(1221 , it has been recently argued [5j that one should com- 
pare G(w) with l/(e(w) — n 2 ), the so-called "rotational 
modulus" . This is true if the equation is valid; it is not 
true if one has a different decay mechanism of the dielec- 
tric polarization. Take, for instance, the Debye model of 
eq. (|17|) . The true relaxation peak lies at 1/td in c(lo) as 
it should, but shifts with varying Ae V ib in the rotational 
modulus l/(e(w) — n 2 ). 

An alternative to the Debye-Gemant-DiMarzio-Bishop 
scheme is a recent proposal |7J for the structural relax- 
ation in the energy landscape |20l (for a review of the 
energy landscape concept se e \2a ) The proposal is based 
on the experimental finding [13 of two time scales in the 
undercooled liquid. The faster one is the shear stress de- 
cay, described by g(u). The slower one is the decay of 
the structural potential energy, seen in dynamic heat ca- 
pacity measurements or in transient grating experiments 

To understand this behavior, let us assume that the 
shear stress decay occurs via thermally activated pro- 
cesses in the energy landscape |20j, |2l|, without speci- 
fying the exact nature of the thermally activated pro- 
cesses. They might be "flow by shoving" [2^] or some 
other mechanism [23|; for our purpose, it suffices to as- 
sume that the shear decay happens by repeated passages 
of the system from one minimum of the energy landscape 
to another. The energy landscape picture provides a nat- 
ural explanation for the dynamical heterogeneity seen in 
numerous experiments |24| . One can show that an initial 
and a retarded part of the flow process follow from the 
energy landscape concept under reasonable assumptions. 

To see this, consider two structural rearrangements oc- 
curring in different parts of the sample. Each of them can 
be characterized by an asymmetric double-well potential 
in the corresponding configurational coordinate, leading 
from one potential energy minimum to the other. Since 
the two structural rearrangements are far apart on a mi- 
croscopic scale, the two minima of the first structural re- 
arrangement change only slightly by a jump in the other 



G(t) 
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FIG. 2: Spring model for the decay of the structural potential 
energy of an undercooled liquid. 
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double- well. 

In a macroscopic sample, there is a very large number 
of possible structural rearrangements. Each of them is 
a local process, and each of them influences every other 
one. 

The consideration shows that one has two relaxation 
times in each double- well potential, a lifetime of the pop- 
ulation of the two minima and a lifetime of the double- 
well potential itself. For the latter one, the Maxwell 
time is crucial, because it sets a lower boundary for the 
time scale on which the core region of the relaxation can 
change its shape. If the core region can change its shape, 
the core relaxation can disappear or make way for a dif- 
ferent relaxation. This relaxation life time is not deter- 
mined by a single jump, but by the joint effect of all 
the other relaxations in the sample, a diffusional motion 
through phase space. The dielectric data seem to see 
this final relaxation time. In fact, the dielectric relax- 
ation time in glycerol and OTP is close to the spin-echo 
relaxation time at the first sharp diffraction peak , the 
decay time of the short range order. 

Naturally, the energy landscape itself does not change 
in time; nevertheless, we believe that the concept of a 
decay of a local double-well potential, for simplicity as- 
sumed to be the same for each double- well in the system, 
catches an essential feature of the motion of the system 
in the fixed energy landscape. Qualitatively, such a pic- 
ture is also compatible with NMR findings 25] in toluene 
and glycerol, showing a small number of larger-angle ro- 
tational jumps (attributable to the hopping part) and 
a larger number of small-angle rotational jumps of the 
molecules (attributable to the diffusional part). 

The final part of the equilibration is not (or maybe only 
partially) seen in the stress decay, because it can only 
happen after the stress decay. As long as there is still 
some appreciable shear rigidity, higher and higher barri- 
ers are jumped over in order to equilibrate the remain- 
ing average stress to zero. This mechanism determines 
the Maxwell time, and a corresponding Maxwell barrier. 
Roughly speaking, the part of the motion through phase 
space due to the initial change of conditions changes from 
hopping to diffusion after the Maxwell time. If an en- 
ergy barrier is higher than the Maxwell barrier, it is not 
jumped over, but it flows away. 

Here, we follow ref. in assuming that the retarded 
part of the process can be described as the decay of an 
energy stored in a harmonic spring rG<x> in series with a 
time-dependent spring G(t) (see Fig. 2). The lower r, the 
more retarded is this final stage of the equilibration, r is 
the ratio of the energy stored in G(t) to the one stored in 
the harmonic spring, which does only decay by flowing 
over into the spring G{t). 

Note that not the structural potential energy itself 
decays, but its difference to an average temperature- 
dependent value determined by the entropy of the inher- 
ent structure of the energy landscape. Our assumptions 
imply that the difference between the potential struc- 
tural energy of a given minimum of the energy landscape 



and the average value has always the ratio r between its 
harmonic part and its long-range stress component. 

In the decay of the structural potential energy of the 
undercooled liquid, one expects to see both the initial 
hopping and the final diffusion. The previous paper |jj 
took only this second retarded action into account. But 
naturally one has to take the initial process into account 
as well. Each jump in an asymmetric double- well changes 
the structural potential energy. This forces one to intro- 
duce an additional dimensionless parameter /, the frac- 
tion of the structural potential energy which equilibrates 
via the initial hopping. In a globally connected energy 
landscape pfj. this parameter / would be one, but in a 
real energy landscape it is expected to be considerably 
smaller. 

We assume that this initial part has the same time de- 
pendence as the shear modulus G(t). Then the Fourier 
transform of the normalized decay function of the struc- 
tural potential energy is 



$(w) = /flM + (l-/): 



1 +r 



75H ' 



(23) 



where g(u>) = G(w)/G os . 

If the dielectric polarization has indeed the same time 
dependence as the structural potential energy difference, 
we have to insert this $(w) into eq. Ijl2|l and obtain 



5(0)- 



+ (1-/) 



1 - gH 

l + fl (w)/r" 



(24) 



Experimentally, one finds the loss peak of the di- 
electric constant close to the one of the heat capacity 
El EiJ H3 In fact, one can argue that the elec- 
tric field introduces an imbalance between the different 
structural realizations of the undercooled liquid at the 
given temperature. If one mirrors the sample at a plane 
perpendicular to the electric field, one gets a state with 
the same structural potential energy, but with opposite 
electric dipole moment. These two states have the same 
energy in the absence of the electric field, but a different 
energy in the electric field. Thus the switch-off of the 
field leaves an imbalance in the structural entropy. It 
should equilibrate by the same mechanism as the imbal- 
ance introduced by a temperature jump. On the basis 
of this argument, one understands the coincidence of the 
dielectric relaxation time with the dynamic heat capacity 
and the transient grating one. 

If the decay mechanism of the dielectric polarization 
is not the diffusion of a single molecule, but the pas- 
sage from one of the minima of the energy landscape to 
another, one needs no longer distinguish between the in- 
ternal field at the molecule and the externally applied 
field, because the externally applied field acts directly on 
the electric dipole moment difference of the two minima; 
all interaction effects are already included in this dipole 
moment difference. The complicated many-b ody prob- 
lem of a system of interacting electric dipoles [T(| enters 
only via the temperature dependence of e(0). 
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It is interesting to consider the limiting cases of very 
small and very large r in the retarded term of eq. (|24p. 
For very small r, the denominator dominates the behav- 
ior and we approximate the Gemant-DiMarzio-Bishop 
case of eq. (|22|l . For very large r, the harmonic spring is 
like a rigid connection and the relaxation peaks of G(w) 
and e(w) become identical. Note that this limit can never 
be attained by the Gemant-DiMarzio-Bishop extension of 
the Debye scheme. There, it corresponds to the limit of 
very small c r , in which the relaxation peak disappears, 
leaving only the vibrational component. 



III. DIELECTRIC, THERMAL AND 
ROTATIONAL RELAXATION TIMES IN SALOL 

For the comparison of dielectric and thermal relaxation 
times, we add a new example, salol (glass temperature 
213 to 217 K, mode coupling critical temperature 255 to 
265 K), to the three examples glycerol, propylene car- 
bonate and OTP discussed in the previous work Q. All 
measured relaxation times are compared to the Maxwell 
time, calculated from fits of the infinite frequency shear 
modulus and of the viscosity. The infinite frequency 
shear modulus is taken from light scattering Brillouin 
data. As it turns out, it is better to use longitudinal 
sound wave data, relating the transverse sound velocity 
v t to the longitudinal vi by vi/v t ~ 1.8. The mistake of 
this approximation is small on the scale of factors of ten 
considered here. The advantages are (i) the longitudi- 
nal infinite frequency sound velocity is much more easily 
determined (ii) one has data up to the highest tempera- 
tures. To take an example, G^ from transverse Brillouin 
OTP data (2(| extrapolates to zero at 348 K (not at 308 
K as stated erroneously in the previous paper), while one 
still has data to compare up to 380 K. 

In salol, the density follows the relation (2?J 



p= 1451.6 - 0.857T 



(25) 



with p in kg/m 3 and T in Kelvin. The longitudinal infi- 
nite frequency sound velocity |28j | 



(26) 



with the glass transition temperature T g — 21%K. With 
our recipe for the ratio vi/v t , this yields 



Goo — 2.2 ( -p^ 



GPa. 



(27) 



The viscosity (2^, y(J, |31| is parametrized in terms of 
two overlapping Vogel-Fulcher laws 



log r] = logr/oi + 



B, 



T - Toi 



(28) 



with i = 1 and i = 2, respectively. The first of these 
two is valid below a temperature T\ , the second above a 



temperature T 2 < T±. Between T 2 and T\, one takes a 
linear interpolation between the two to ensure continuity. 
For salol, log 7701 = —36.9 and log 7702 = —3.68 (77 in Pas), 
Bi = 7500 K and B 2 = 151 K, T 01 = 59 K and T 02 = 225 
K, T x = 256 K and T 2 = 248 K. 

Figure 3 compares dielectric j^E^, nea * capacity |32^ 
and transient grating [35l | relaxation times to rotational 
relaxation times extracted from longitudinal Brillouin 
light scattering data j34[. The times were either recal- 
culated (3(| or refitted in terms of a Kohlrausch function 
expi-it/TKww) 13 ). 

Fig. 3 contains no mechanical shear relaxation 
times, but one knows that these have to lie below the 
Maxwell time, because for a shear modulus following the 
Kohlrausch function exp(— (t/ T shear) ) 



T~she 



T~Ma 



(29) 



Usually, (3 lies between 0.4 and 0.6, so the ratio should be 
between one third and two thirds; the shear Kohlrausch 
relaxation time should be a factor 1.5 to 3 shorter than 
the Maxwell time. 

Fig. 3 shows once again that the rotational times keep 
close to the entropy relaxation times in the heat capac- 
ity and in the transient grating measurements, thus sup- 
porting the proposed energy landscape mechanism for 
the relaxation of the molecular orientation. In principle, 
the transient grating measurements should also be evalu- 
ated in the Ansatz of Pick and Dreyfus pV?! l3^ | , because 
the signal comes from a grating of both temperature and 
local molecular orientation. However, an experimental 
separation 39] of the contributions demonstrated again 
the equality of the two relaxation functions within exper- 
imental error. 

Note that the method developed by Dreyfus and Pick 
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FIG. 3: Kohlrausch- Williams- Watts relaxation times in sa- 
lol, normalized to the Maxwell time as described in the text. 
Symbols: pluses dielectric |3^ |: crosses dielectric pisl ]: full tri- 
angles heat capacity |3^|: full diamonds rotational relaxation 
times extracted from Brillouin light scattering |34| ; full circles 
transient grating data |3q . 
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[371 l38| to determine rotational relaxation times from 
light scattering Brillouin data seems to work remarkably 
well. According to their Ansatz, one should see the sec- 
ond Legendre polynomial. For continuous rotational dif- 
fusion, their relaxation times should be a factor of three 
smaller than the dielectric ones. One does not see this 
factor of three; the two sets of data are remarkably close 
to each other over ten decades in relaxation time, even 
making the same wiggles compared to the Maxwell time. 
In the energy landscape mechanism, one does not expect 
the factor three; there, it depends on the jump angle dis- 
tribution which can be different for different glass formers 

Hi. 

The marked temperature dependence of the ratio be- 
tween rotational relaxation time and Maxwell time indi- 
cates a temperature dependence of the retardation pa- 
rameter r in eq. (|24|) . The rise at low temperature is 
also seen in glycerol and in propylene carbonate Q ■ The 
salol kinks at 260 and 240 K have been also seen in the 
"Stickel plot" |3i| of the dielectric data (a plot of the 
inverse square root of the temperature derivative of the 
peak frequency versus temperature, where one sees the 
deviations from a perfect Vogel-Fulcher behavior). 

From the point of view of the mode coupling theory 
[4fH |. the kink at 260 K should mark the crossover from 



an energy landscape behavior below T c to simple liquid 
behavior above. But Fig. 3 shows that the separation 
of time scales persists into the liquid domain in salol. 
This time scale splitting is not the two-stage scenario of 
the mode-coupling theory, because both time scales move 
together with the Maxwell time. In fact, in ref. 0] 
this time scale splitting was also observed for propylene 
carbonate and discussed in terms of the mode coupling 
theory. The a-process of the theory was not attributed 
to the slower, but to the faster process. 



IV. SUMMARY 

The textbook relation between dielectric polarization 
decay and dielectric constant is applied to the simple 
Debye case, to its extension to viscoelasticity and to a 
newly proposed energy landscape mechanism. 

Data in salol show once again the equality of dielectric 
and structural relaxation times, supporting the energy 
landscape mechanism. 

We thank Catherine Dreyfus and Robert Pick for help- 
ful discussions. 
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